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ABSTRACT 



The application is reported of a United Sensor Corpor- 
ation 5 -hole pressure probe to determine the flow field 
upstream and downstream of the rotor of a transonic com- 
pressor operating to near sonic speeds. Calibration of the 
probe was reduced to a set of equations to allow off-line 
data reduction and presentation .using a Hewlett-Packard 
Model 983OA Calculator system. The flow velocities and 
angles measured upstream of the rotor were compared with 
results of compressible flow calculations. Blockage factors 
were determined which related analytical approximations 
of the flow incident at the rotor to the measured flow 
rate. The measurements downstream of the rotor were com- 
pared with those of a second probe of different design. 

Both probes required corrections to their calibrations 
which were established in an axial annulat* flow of sim- 
ilar geometric cross section. It is shown that the results 
derived from the corrected calibration of the 5 -hole probe 
were the more accurate. 
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I. INTRODUCTION 



The study reported here is an initial part of a program 
to obtain a full understanding of the flow within a tran- 
sonic compressor stage, with the ultimate object of obtaining 
improved designs. Toward this goal, a transonic compressor 
of the "impulse" type was designed and built by Dr. M.H. 

Vavra in the Turbopropulsion Laboratories at the Naval 
Postgraduate School. The work v/as sponsored by Naval Air 
Systems Command, Code 310 through the office of Mr. H.J. 
Mueller. The so-called TRANSX compressor has been operated 
to of design speed. 

The inlet of the TRANSX compressor is designed to con- 
tract and to rapidly accelerate the flow in order to -gener- 
ate near uniform conditions, with thin boundary layers, at 
the rotor face. It is necessary to determine the radial dis- 
tributions of the flow velocity and angle in order to estab- 
lish the operating condition and to evaluate the performance 
of the rotor blading. At or near design speeds, the possibil- 
ity that wakes in the flow might induce vibrations in the 
highly stressed rotor blades, prohibits the use of probes 
just upstream of the rotor. Ttierefore it is necessary to be 
able to predict velocity and pitch angle distributions at 
high speeds from data gathered at lower speeds. Reference 1 
gives analytical approximations of the rotor face flow velo- 
city and angle distributions which were obtained from com- 
pressible flow calculations. The first goal of the present 
study was to determine the flow into the rotor by measurement, 
and to compare the results with the analytical predictions. 
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A United Sensor 5-hole cylindrical pressure probe was selec- 
ted for its ability to measure pitch angle and for ease of 
mounting. The probe calibration was first established in a 
uniform jet. Following a description of the TRANSX compress- 
or in Section II, the investigation upstream of the rotor is 
described in Section III of this report. 

Measurement of the flow downstream of the rotor is 
necessary in order to establish .blade element performance. 
The calibrated United Sensor 5-hole probe and a NASA combin- 
ation probe [Ref, 2] were used simultaneously to conduct 
radial surveys of the flow downstream of the rotor. In this 
manner redundant measurements of Mach number and flow angle 
were obtained using two completely independent probes and 
their respective calibrations. It was found that the Mach 
numbers evaluated by the two probes disagreed by as much as 
10^ when the calibrations, established in the uniform jet, 
were applied to the probe measurements made downstream of 
the rotor. Methods were devised for correcting the calibra- 
tions established in the uniform jet for application of the 
probes in the smaller compressor annulus. Following the 
application of corrections, which were greatest for the NASA 
probe, the disagreement in Mach number profiles was greatly 
reduced; in some cases to less than Vfo, The measurements 
carried out downstream of the rotor using the 5-hole probe 
are described in Section IV, 

The calibration of the 5-hole probe and the method de- 
vised for correcting the calibration for application in the 
compressor annulus are described fully in Appendix A, 
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Inter- re la ted non-linear calibration relations are therein 
represented by analytic functions which allow computer re- 
duction of probe survey measurements. The methods described 
in Appendix A can be adapted to other probes and to measure 
ments with the same probe in passages of different geometry 

II. THE TRANSONIC COMPRESSOR 

A. DESCRIPTION 

A schematic of the compressor test rig is shown in Fig. 
1. The test rig consists of an air- turbine drive unit and 
an induction section containing a filter, a throttle, a 
settling chamber and a flow measuring nozzle. The air tur- 
bine is designed to deliver 450 H.P. at 30,000 RPM. The 
transonic compressor is designed to operate at 30,460 RPM 
with a relative tip Mach number of 1.5» The compressor flow 
rate is controlled by a rotating throttle plate T]« 

The test rig, located in an explosion proof cell, is con- 
trolled by an operator at a remote console, A complete des- 
cription of the laboratory facilities sre given in Ref, 3» 
and the details- of the test rig design is given in Ref. 4, 
Figure 2 shows the compressor assembly and identifies the 
instrumentation stations 0 through 3 referred to in the 
following discussion. As indicated in view J of Fig, 2, 
peripheral surveys were made by .rotating the outer casing, 

B. INSTRUMENTATION 

The transonic compressor is equipped with stationary 
probes to provide routine performance data. The measurement 
of torque was not available at any time during this 
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investigation, since the honeycomb flow straightener (Fig, 

3) was removed to allow higher flov/ rates. The fixed instru- 
mentation is not discussed in this report. 

The results presented here were obtained using a United 
Sensor DA-125 three dimensional pressure probe. The probe 
calibration is reported in detail in Appendix A, The probe 
mounting arrangement is shown in Fig, 3* Radial surveys were 
•taken at a fixed peripheral angle by manually positioning 
the probe using a traverse unit. Peripheral surveys were 
made by manually rotating the outer casing, causing the 
probe to move through an arc at a fixed radius. All press- 
ures were measured by a calibrated transducer in a scani- 
valve arrangement. Predetermined angular and radial posi- 
tions were input into the data logging system on thumb 
wheels, 

C, DATA ACQUISITION AND REDUCTION 

A data logging system [Ref, 3] was used to scan and re- 
cord pressure, temperature and position data. The recorded 
output was on punched paper tape. Reduction of data was 
carried out using a Hewlett-Packard Model 98JOA Calculator 
system, shown schematically in Fig, 4, 

A program, designated TCR Reduction (Station 1), was 
written to reduce data from rotor inlet flow measurements. 
Data were input from the punched tapes. The program calcula- 
ted and printed the non-dimensional radial probe position, 
compressor case angle, flow pitch and yaw angles, non-dimen- 
sional velocity, velocity and Mach number. These results 
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were stored on digital magnetic tape and the radial distri- 
butions were graphed using the X-Y plotter (Fig, 4), 

A program entitled TCR Reduction (Station 2) was 
written to calculate the flow downstream of the rotor from 
probe measurements. The program accepted data from magnetic 
tape, rather than punched paper tape. The paper tape data 
were therefore transferred to magnetic tape before beginning 
data reduction. This procedure was selected because data 
from two different probes were required to be in the same 
form to allow distribution plots to be made on a common 
graph, V/hen data had been reduced, the calculated results 
were stored on magnetic tape and the radial profiles were 
produced on the X-Y plotter, 

III, FLOW INTO THE ROTOR 
A. DESCRIPTION OF TESTS 

The purpose of these surveys was to measure the flov/ 
properties at the rotor inlet and to compare the* results 
with analytical distributions reported in Ref. 1. Prelimin- 
ary impact pressure surveys were made prior to this investi- 
gation, and as a result, a bellmouth inlet fairing was added 
(Fig, 2) to improve inlet profiles. 

The 5-hole probe was positioned ahead of the rotor at 
station 1, The flow profiles were measured in two separate 
runs (numbered 33 and 39) • All surveys were radial, from 
hub-to-tip. The speed was set at 50?^ of design speed in run 
33* and varied from 65?^ to 4o^ of design speed in rUn 39* 
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Throttle settings corresponded to fully open, near surge, or 
design point for each speed of a given survey. The case 
angle, or peripheral location of the probe, was not consider- 
ed to be significant for measurements at this station. 

The surveys began with the sensing holes 0,2? inch 
from the rotor hub. The probe was then withdrawn in 0,1 
inch intervals until the total pressure sensing port was 
flush with the inner surface of the casing. At each point, 
pressure, temperature and position data were recorded on 
punched tape. These data were reduced to velocities, angles 
and non-dimensional position, 

B, RESULTS 

Radial distributions of the pitch angle and Mach number 
derived from probe measurements are shown in Fig, 5 to Fig, 

11, The measured data were reduced using the probe calibra- 
tion described in Appendix A, section A2, The radial distri- 
bution of pitch angle, derived in Ref, 1, is also shown on 
the plots. The data were taken in two separate tests. Re- 
sults from run 33 in Figs, 5-7 show good agreement with the 
derived profile except at the upper and lower 20^ of the dis- 
tribution, Note the wall effect on pitch angles for radial 
positions greater than R=0,85* Within the measurement accur- 
acy, distribution of pitch angle is independent of flow rate 
for all surveys taken in a given test. The results of run 39 
(Figs, 9-12) show an overall reduced level of pitch angle, but 
also shows no dependence upon throttle position or speed. 

The reason for the reduced levels is not known and should 
first be confirmed by further measurements. 
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In order to plot the analytical Mach number distribu- 
tion as a comparison with the measurements plotted in Fig, 5 
to Fig. 11, the area blockage factor must first be determin- 
ed, Therefore the measured distribution of the mass flux com- 
pared with the analytical distribution was chosen to repre- 
sent the same information, and is presented in the following 
section. 

The calibration correction- for application in the com- 
pressor annulus does not apply to measurements at the rotor 
inlet. However, this correction was applied to the profile 
in Fig; 6, and the results displayed in Fig, 8, The correc- 
tion applied near the outer wall causes the profile to 
better match the analytical distribution at the wall. Since 
the wall angle at the inlet is the same as that in the cali- 
bration annulus, the effect brought about by applying the 
correction supports its validity for this region, 

C. DETERMINATION OF BLOCKAGE FACTOR 



Details of blockage factor calculations are given in 

Appendix B, Reference 1 determined a viscous blockage factor, 

Kp =0,989* which accounted for wall boundary layer displace- 
^1 

ment thickness. Figures I 3 and l4 shov/ distributions of the 
non-dimensional mass flux, or flow function obtained from 
probe measurements and a comparison with the analytical dis- 
tribution given in Ref, 1 and Appendix B, It is shown in 
Appendix B that by equating the integral of the measured mass 
flux, over the region for which the probe measurements are 
free of wall effects, to the integral of the analytical mass 
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flux distribution, a blockage factor can be defined. The 
blockage factors listed in Figures 13 and l4 were obtained 
in this manner and the analytical distributions are consis- 
tent with these values. 

The prediction that the mass flux increases with in- 
creasing radius is supported by the measurements, but the mea- 
sured variation is larger than anticipated. This could be 
the result of the inlet bellmouth velocity profiles which has 
yet to be measured. 

The mass flux profiles are quite similar in shape, but 
considerable variation occurs in the blockage factors de- 
duced from the measurements (,963 to .984), This is possibly 
due to the measurement of the total mass flow rate, which 
was observed to vary during the surveys. As shown in Appen- 
dix B, an attempt was made to eliminate this effect; however 
insufficient data was taken to make the process statisti- 
cally accurate. It is observed that the smallest blockage 
factors were obtained at the lowest flow rates, for which 
the measurements are least accurate. However, it is' also 
noted that at highly reduced throttle settings, greater dis- 
tortion appears in the profiles near the rotor tip, suggest- 
ing that an increase in blockage does indeed occur. This 
question may be resolved by collecting more data witli the 
objective of determining how blockage varies with RPM and 
throttling conditions. 

Figures 15 and I 6 shov; calculated flow distributions 
into the rotor using the method of Ref, 1 with a blockage 
factor of 0 . 98 , Note, in Fig. 15i that choking of the inlet 
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occurs at slightly greater than 21 Ib/sec, As blockage in- 
creases (decreasing effective cross sectional inlet area) 
choking occurs at lower flow rates. Figure l6 shows that the 
rotor blade relative Mach number is subsonic everywhere at 
the speeds encountered thus far. At design speed (30»460 
RPM) 80 ^ of the blading will see supersonic relative veloci- 
ties, assuming a blockage factor of approximately 0,98, 

I 

IV, FLOW DOWNSTREAM OF THE ROTOR 
A, DESCRIPTION OF TESTS 

The United Sensor probe and a NASA combination probe 
[Ref, 2] were mounted at station 2, separated peripherally 
"by a 45 ° angle. There was a 1,5° difference in the position 
of each probe with respect to a stator blade leading edge. 
Peripheral surveys were plotted such that peripheral angle 
displacement of each probe was measured with a stator blade 
leading edge as reference. During radial surveys, the two 
probes were moved outwards and data from both were taken 
simultaneously. 

Three separate runs, numbered 43, 44 and 45 are re- 
ported, The tests consisted of one radial survey and one 
peripheral survey at the passage centerline for each of the 
three runs. Runs 43 and 44 were conducted at the same RPM, 
but .at different throttle settings. Runs 44 and 45 were con- 
ducted at the same throttle setting, but at different 
speeds. 

The purpose of these surveys was to verify measurements 
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through redundancy and to calculate rotor performance. 
Agreement was obtained only after probe calibrations were 
corrected, as described in Appendix A, The verification of 
the measurements themselves is described in this section, and 
the procedure for calibration, correction and application to 
the 5-hole pressure probe is described in Appendix A. The 
evaluation of rotor blade element performance will be re- 
ported separately. 

B, RESULTS 

Figure 17 shov/s a comparison of Mach numbers for run 45, 
point 3» evaluated from 5-hole and NASA probe measurements 
using calibrations established in a uniform jet. Approxi- 
mately 10 ^ disagreement between the measured Mach numbers 
was observed. Since neither the static pressure distribu- 
tions calculated from the 5-hole probe calibration nor that 
from the NASA probe calibration were consistent with the 
measured compressor hub and tip wall pressures, the Mach 
number distributions shown in Fig, 17 were not acceptable. 

It had also been observed from surveys of the uniform jet, 
described in Appendix A, that the probes indicated uniform 
flow profiles over different regions of the flov/. In parti- 
cular, the NASA probe failed to indicate uniform conditions 
until it was traversed radially more than 3 inches into the 
flow. The effect of "immersion" or "wall effect" on the 
calibration of the two probes was the probable cause of the 
observed disagreement exhibited in Fig. 17, Since the magni- 
tude of the immersion effect, in the uniform jet, was quite 
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different for each probe, the approach taken was to first 
devise separate corrections to the calibrations of the two 
probes. Agreement between the corrected survey measurements 
with both probes, and consistency of the calculated static 
pressure distributions with the measurements at the hub and 
tip walls v/ere sought. The corrections v/ere based upon sur- 
vey measurements made in an annular axial flow passage which 
• was similar in size to the compressor annulus downstream of 
the rotor. The procedure for correcting the 5-hole probe 
calibration is described in Section A3 of Appendix A, 

The results shown in Fig, 18 to Fig, 2° were obtained 
follov/ing the application of corrections to the measurements 
of both probes. The upper sections of Figs, 18-23 show the 
impact pressure (to the right) and indicated static press- 
ure (dashed lines) measured by the probes. The solid lines 
tov/ard the left are the distributions of static pressure 
evaluated from the probe measurements after the annulus 
corrections v/ere applied. The lower sections of Fig, 18 to 
Fig, 23 compare Mach number and flow angles derived from the 
probe measurements. Figures 24 to 29 show pitch angle dis- 
tributions from the 5-hole probe measurements. 

The results shown in Figs, 18 to 29 were from the three 
tests mentioned above. Figures 18, 19 and 20 show radial 
surveys at ^ 0 % of design speed and design incidence (deter- 
mined from Fig, 15) • 50?5 speed near surge and 65 ^ speed near 
design incidence, respectively. Figures 21, 22 and 23 show 
peripheral surveys of the centerline of the channel at the 
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same test conditions. Figures 24 to 29 show the resulting 
pitch angle distributions corresponding to tests shown in 
Figs, 18 to 23» 

C,' DISCUSSION 

1, Radial Surveys 

The radial surveys are shown in Figs, 18 to 20, The 
Mach number comparison between the two probes is in much 
better agreement than the uncorrected results of Fig, 17. 

The differences close to the hub are a result of the static 
pressure, rather than total pressure measurements. It was 
noted that corrected static pressure measurements for the 
5-hole probe tended toward measured hub and tip wall press- 
ures, whereas the NASA probe corrected static pressure tend- 
ed to values that were too high. As an example, for the test 
shown in Fig, 20, the measured hub and tip wall pressures 
were -5*3 and 23,1 inches of water respectively. Therefore, 
the Mach number calculated from the 5-bole probe measurements 
is considered to be more accurate than that given by the NASA 
probe , 

Yaw angle measurements of flow behind the rotor are in 
agreement to better than one degree. The reason for the resi- 
dual disagreement is not clear. It was noticed that the 
largest and most consistent disagreement occurred at the 
highest speed and flow rate (Fig, 20), In Fig, 21 the flow 
angle at the rotor tip is large; corresponding to a low axial 
velocity component. This suggests a stalled condition at the 

rotor tips. 

Total pressure measurements do not agree as well as in 
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the steady flow in which the probes were calibrated (Appendix 
A). Leaks were thought to have been eliminated from the 
scanivalve and pressure transducer data system. There is 
a possibility that the total pressure differences were 
caused by the unsteady nature of the flow. This argument is 
supported by Fig. 19# corresponding to a near stall condi- 
tion, for which additional unsteadiness might be expected. 
However, at the higher frequencies associated with the 
higher RPM in Fig, 20, the disagreement is greater. Such 
would not be the case if the increased differences were 
caused by the unsteady response of the probes, 

2, Peripheral Surveys 

Figures 21 to 23 show the results of peripheral 
surveys taken at the passage centerline, corresponding to 
test conditions shown in Figs, 18 to 20, In Figs, 21 and 23, 
the impact pressure difference is small and uniform with re- 
spect to angular position. (For these tests, profiles 
hub-to-tip, were also similar), At the near stall condition 
in Fig, 22, a varia-fen in impact pressure v^ith peripheral 
position is evident. The two measurements agree only over a 
small portion of a blade space. 

The peripheral variation in corrected static pressure, 
as indicated by both probes, was quite pronounced. With the 
stagnation pressures being relatively constant, a variation 
in Mach number proportional to the probe static pressure 
differences v/as measured. Inspection of these figures leads 
to the conclusion that the stator bladesvere affecting the 
flow upstream. The minimum Mach number and flow angle 
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(Fig. 22) was measured in line with the leading edge of a 
stator blade, at a case angle of approximately 363°» That 
angle is consistent with the geometry of the flow into the 
stator. It is clear that peripheral variations in flow pro- 
perties must be taken into account for specification of 
rotor performance from probe measurements at station 2, 

This is particularly true at reduced throttle settings, as 
in Fig. 22. 

The distribution of static pressure from the NASA 
probe shown in Fig, 22 is similar in shape to those from both 
probes in Figs, 21 and 23, The more uniform distribution of 
the 5-hole probe static pressure is a marked departure from 
the distributions of Figs, 21 and 23, Since the probes were 
separated by 3 stator blade passages, it is possible that the 
stator passage downstream of the 5-hole probe contained a 
separation, or local stalled condition, whereas the passage 
downstream of the NASA probe was unaffected. 

3. Pitch Angle Measurement 

The pitch angle distributions measured downstream of 
the rotor are shown in Figs, 24 to 29 and correspond to the 
surveys of Figs, 18 to 23, The variation of pitch angle with 
radius v/as nearly linear at design throttle settings (Figs, 

24 and 26), The pitch angle is indicated to be approximately 
-1° close to the outer wall, as might be expected due to 
boundary layer effects. These measurements, which rely upon 
the annulus corrections described in Appendix A, are remark- 
ably consistent with the expected distribution. At the point 
where probe measurements were taken, the hub wall angle was 
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+12°, Extrapolation of the angle distribution to the wall, 
in Figs, 24 and 26, is consistent with this value. The re- 
sults in Fig, 25 , for the near surge condition, shows a re- 
gion of nearly constant pitch angle near the inner wall and 
a tendency to increase rapidly toward the 12° hub wall angle 
This is consistent with rake measurements of impact pressure 
taken at the stator exit, which indicated a region of sep- 
arated flow at the hub wall, Th-is would have the effect of 
creating a large pitch angle gradient near the hub. The lar- 
ger region of negative pitch angle near the tip would be con 
sis tent with an increase in wall boundary layer thickness. 
This region also exhibits a large pitch angle gradient when 
compared to the center of the distribution. 

Peripheral variations in pitch angle, shown in Figs, 

27 to 29 indicate little change at or near design cond- 
itions, However, near surge (Fig, 28) the flow angle incr- 
eased on one side of the stator blade and decreased on the 
other side. Though it appears to be a significant effect, 
the variation is no more than one degree either side of 
the mean line and v/ould not be significant in rotor per- 
formance calculations. 
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V. CONCLUSION 



. In this study a 5-hole pressure probe was used success- 
fully to determine the distribution of velocity and flow 
angle upstream of a compressor rotor. Blockage factors were 
determined from measurements, which allow the use of analyti- 
cal expressions to relate the measured flow rate to the in- 
cident flovz properties at the rotor. 

The same probe v/as applied downstream of the rotor. 

The measurements v/ere compared to those of a NASA probe, 
which did not provide a measurement of pitch angle. Part of 
the differences in the Mach number distributions, seen near 
the hub, was probably for this reason. The 5-hole probe 
gave more accurate Mach number distributions mainly because 
the static pressure corrections for passage effects were, 
smaller than those for the NASA probe. 

The differences in the radial distributions of impact 
pressure were partially the effect of the difference in the 
location of the probes with respect to the nearest stator 
blade leading edge. Although that difference was only 1.5^* 
at decreased throttle settings the peripheral variations in 
the indicated flow profiles became significant. In future 
tests radial surveys should be conducted with both probes in 
the flow simultaneously, at fixed speed and throttle setting; 
then the probes should be interchanged and the surveys re- 
peated without changing conditions. The resulting plots would 
differ only as they are affected by the slight peripheral 
displacement of the probes. Peripheral survey variations must 
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be accounted for if probe measurements at station 2 are to 
be used for derivation of rotor performance. 

The effects of unsteady response on probe measurements 
should also be examined. It is reported that the average in- 
dicated total pressure in unsteady flow can be greater than 
the time average of the total pressure [Ref. 6]. This implies 
that two different probes in the same steady flow could 
measure identical total pressures, and when placed in un- 
steady flow might disagree as a function of speed and throt- 
tle setting. While using Reference 6, this effect was calcu- 
lated to be small for the probes and conditions in the present 
study. The differences in the measurements suggest that this 
conclusion should be re-examined. 

The correction devised to compensate for wall effect on 
the 5-hole probe pitch angle measurements v/as very effective 
judging from the plotted pitch angle distributions at station 
2. The overall effect of correcting the 5-hole probe calibra- 
tion for application to the specific flow geometry of the 
transonic compressor, points out the necessity of Careful 
calibration of pressure probes used for making measurements 
of internal flov/s. 
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Figure 1, The Compressor Test Rig 
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Figure 2, The Compressor Assembly 
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Figure 4 , Schematic of the Hewlett-Packard 
Model 983OA Calculator System 
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FIGURE S. PITCH RUBLE flUl> MHCH NO. DISTRIBUTION RT 5T.qTinN 
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FISURC B. PITCH HHELE RN1> MHCH ND. IM5TRIBUTIDN RT STHTION 
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FIGURE 7. PITCH RNGLE RND MRCH ND. I> 1 5TR I HUT I DN HT STRTIDN 
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FIGURE B. PITCH RNGLE RND MHCH NO. DIGTRISUTIDN RT STBTIDN 
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F1BURE:3. pitch angle and MHCH Nn. I>I5TRIBUTinN RT STRTICN 1. 
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figure: 10. pitch angle RND MRCH ND. DISTRIEUTION RT STHTIRN 1. 
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FIGURE I I. PITCH RNHLE RND MRCH ND, DISTRIBUTION RT 5TRT10N 
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FIGURE 12. PITCH ANGLE AND MACH ND. l>l5TRIBUTinN AT STRTIDN !. 
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riE.t3. D15TRIBUT1DN DF FLOW FUNCTION RHFflD OF THE ROTOR. 
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FIG. IM. DISTRISUTIDN OF FLDH FUNCTION RHEflI> OF THE ROTOR. 
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FIG. 15. TRRNSONIC COMPRESSOR - CRLCULRTED FLOW INTO THE ROTOR 



COMPRESSOR RPM 




CK 



H4 



a 

lij 

to 

:c 

Z 

X 

E 

Z 

«a uj 
j > 

H- 

0; 

Ui 

E 



Lrl 

ts 



ta 

IS 



47 



FIE. IE. TRHMSDNIC COMPRESSOR - CRLCULRTED FLOW INTO THE ROTOR 
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Figure 17. Mach Number Distribution Prior To 
Application of the Annulus Corrections 
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FIGURE 2M. PITCH RNGLE I> 1 STR 1 BUT J UN BEHIND THE RDTDR. 
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FIGURE 25. PITCH ANGLE DISTRIBUTION BEHIND THE ROTOR. 
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FIGURE 2H. PITCH BNGLE MSTRIBUTIDN EEHIND THE RDTDR. 
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FI EURE 27. PITCH BNELE DISTRIBUTION BEHIND THE ROTOR 
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r I EURE 20. PITCH BNSLE DlSTRIgUTlDN BEHIND THE RDTDR 
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FlBUf?E2S. PITCH RNSLE DISTRIBUTION BEHIND THE ROTOR. 



Appendix A : Determination of Velocity Vector From 

^"Hole P r essure Probe Measurements 

A .1 Method of Approac h 

The limited axial space available for instrumentation in most turbo- 
machinery flow passages severely limits the types and sizes of pressure 
probes which can be used. The United Sensor DA-125j 5~hole pressure probe 
is designed for this application. The probe, when calibrated, will 
measure flow velocity, yaw angle and pitch angle. A major drawback of 
the probe has been the tedious labor required to extract velocity and an- 
gle data from non-linear, and inter-related calibration curves. The 
purpose of this appendix is to document a method for calibrating the 
United Sensor 5-hole probe such that, in application, data reduction can 
be performed by computer at the time of acquisition. 

Results are given here for a particular probe which was used to 

measure profiles in the TRANSX. compressor. It was found to be necessary 

to correct the measurements made in the compressor annulus for effects 
attributable to the passage size and shape. The method is therefore in 
two parts. First, the probe was calibrated with reference to a Prandtl 
probe in a jet contained within a larger pipe. . Analytical approxima- 
tions were found which were represented as a function of the radial position 

in the annulus. These steps are described in section A2 and section A3 

respectively, where methods of applying the calibration with and without 
the passage correction are also given. 



A. 2 Calibration and Application Free of Wall Effects 
A .2.1 Probe Description and Analysis 

The 5-hole probe is shown in Fig. Al. Characteristics of this and 
similar probes are given in the manufacturer's catalogue (8 )• Briefly, 
the probe is rotated in the flow until the pressure indicated by ports 2 
and 3 (p^ ^ equal. The yaw angle (fv) is then read on a vernier 

scale attached to the top of the probe. (Fig. A2). The average of pressures 
p^ and p^ oP'tained pneumatically through a mixing value) is referred to 

as the "indicated static" pressure of the probe. Port IvPj^) measures "in- 
dicated total" pressure and the difference between the pressures p^ and p^ 
is a measure of the pitch angle. However, the relationship of P 22 > 

(Pi “ Pc) actual static pressure, total pressure and pitch angle must 
4 5 

be established by calibration in a known flow. Indicated static pressure is 

known to be more subject to errors than total pressure, particularly when 

the sensing holes are within about ten probe diameters of a boundary. 

The calibration of the probe is expected to depend on flow Mach number. 

It is more convenient however to non-dimensionalize the velocity using the 

"limiting" velocity, V = /2C T (for a perfect gas), when the stagnation 

X/ P 

temperature is a constant in the flow. The use of limiting velocity is 
discussed in Appendix A of reference 7 . The non-dimensional velocity is 
given by 

X = V/V^ A(l) 

where V = IO 9.62 /T ft/sec. for air when the stagnation temperature, T , 

"C "O’ 

is in^R. As shown in reference 7 > the relationships between the static 
and stagnation gas properties become 
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V 

Y-1 



= (1 - x^) 



V 

Y-1 



The relationship between Mach number and X is given by 



2 2 
X = ^ ^ 

1 + 



A(2) 



A(3) 



The calibration of the 5-hole probe can be represented in the form of 
two functions of Mach number (or X) and pitch angle, 0. 

If p and p are the true stagnation and static pressures, then define 



(X, 







and 



(X, 0) 




A(5) 



to be the required calibration functions. 



probe pressures are taken when p^ = p^ 



where it is understood that the 



The value of X is obtained from true static and stagnation pressures 
using Eq. A(2) : 



X = 




1 - [1 - 




V-1 

Y 



A(6) 



where = atmospheric pressure, and the pressure differences can be 
measured by manometers. 
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To calibrate, the probe must be placed in a flow at a known pitch angle 

and the true static and impact pressures measured by a reference probe. 

From the measured pressures, the value of X is calculated from Eq. A(6) and 

the values of F and F^ from Eq. A(4) and Eq. A(5)» However, in applying 
9 'P 

the probe in an unknown flow, expressions relating only 5 -hole probe 
measurements to the unknowns X and 0 are required. Rearranging Eq. A (4) 
and using Eq. A(2), 



p = (^) Fg (X, 0) [1 -f--} = (^) F^(X, 0) [1 - (1-X^) A(7) 



Define 



F^ (X,0) = F^ (X,0) [1 - (1-X^) 



A(8) 



and, similarly 



FI (X,0) = F. (X, 0) [1 - (1-X^) 



0 



0 



A(9) 



The two functions F^ and F^ can be tabulated as functions of X and 0 using 
the values obtained in the calibration for F and F , at each value of X. 

9 0 

If, in addition, it can be shown that = 1*0 "to within the desired 

accuracy over the range of X and 0 of interest, then, from Eq. A (7) 



F’ (X, 0) = 

B 



pi ~ ^23) 
pi ■ ^A) F; 



A(10) 



and similarly 



Pu - p. 



"0 



A(ll) 



f 
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In using the calibration expressed by Eq. A(10) and Eq. A(ll), it is 

necessary to interpolate in the tables established for and to find 

values of X and 0 for which Eq. A(10) and Eq. A(ll) are simultaneously 

satisfied by the measurements. Alternatively, analytical approximations 

can be derived for F' and Fi which can be solved by iteration. The analy- 

9 0 

tical approximation technique was chosen here. 

A. 2. 2 Calibration Tests in a Jet 

The 5-hole probe was calibrated as shown in Fig. A3. An A.S.M.E. flow 
nozzle , 7 inches in diameter , was flanged between two 36 inch lengths of 
'10 inch diameter steel pipe. (The internal contour of the nozzle was chrome 
plated and polished to a maximum roughness of 32 micro inches). Air was 
supplied to the 10" pipe at controlled flow rates from an industrial com- 
pressor ( 3) via a 4 ft. diameter settling chamber and entrance bellmouth. 

A 12" long bundle of -g-" diameter tubing removed any swirl remaining in the 
flow from the plenum. The flow from the pipe downstream of the flow nozzle 
exhausted into the laboratory. 

A Rrandtl probe was used to obtain the reference measurements of 
static and total pressures. The 5-hole and Prandtl probe tips were posi- 
tioned on opposite sides of the flow centerline at a radius of one inch. 

The measuring station was 2 inches downstream of the nozzle exit face. The 
% 

5-hole probe was mounted through an aluminum block as shown in Fig. a4. 

The block allowed the pitch angle to be set in 5° increments from - 10° to 
+ 10° and in 10° increments from - 40° to + 40° with no change in the tip 
location. All pressure differences were measured using U-tube water mano- 
meters. The stagnation temperature was measured using an iron-constantarn 
thermocouple probe in the flow upstream of the flow nozzle. 
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The flow at the measuring plane was first surveyed using the Prandtl 
probe with the 5-hole probe removed. The result in Fig. A5 showed that the 
velocity (x) was uniform over a 4” diam core, and therefore that the flow 
at the reference probe tip and the 5-hole probe tip were the same to within 
measurement accuracy. 

The procedure for the calibration tests was the following: The flow 
rate from the compressor was set and allowed to stabilize. The 5-hole probe 
was then inserted in turn at each desired value of the pitch angle. Pressure 
differences were recorded by hand from the manometer for the Prandtl and 
.5-hole probes. The thermocouple probe reading and barometric pressure were 
also recorded. 

To cover the range required in the probe application, 0 was varied from 
-5^ to + 4Q0 in 7 steps. Nine values of X were set giving dynamic pressures 
between 9 --d 75 inches of water column. If ± C.025 inches in the accuracy 
of reading each manometer column, the pressure differences are measured to 
± 0.05 inches of water. However, unsteadiness in the air supply over the 
recording time period can be expected to reduce to overall accuracy to ± 0.1 
inches water. 

A. 2. 3 Analysis of the Results 

The recorded data is given in Table Al. The ratio (p^^/p^) was calculated 

to be within 0.2^o of unity over the complete range of X and 0 measured and 

hence was set equal to unity in representing and applying the calibration 

results. Using the data in Table Al, values of F^ and F^ were calculated 

using Eq. A(8) and Eg. A(9)j where F^ and F^ v/ere obtained from Eq. A(4) and 

Eq. A(5) and X was calculated using Eq. A(6). The remaining task was to 

find analytical approximations for the values of F’ and F^ as functions of 

0 0 

X and 0, and to devise an iterative coii^^uter solution. 
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F' and were plotted as functions of 0 for each value of X. The 
results are shown in Fig. a 6 and Fig. A?. Using polynomial regression on 
the Hewlett-Packard 9330A calculator, F' at constant X was satisfactorily 

e 

fit by a Uth order polynomial in 0. F^ at constant X was satisfactorily 
fit by a 2nd order polynomial in 0. This representation however was found 
to be unsuitable in solving iteratively for X and 0 with F^ and F^ given. 
Instead, expressions were found for 



= FjJ (0) 



const 



as in Fig. A7j and 



('■9) l0 = const 



Fig. A8 shows the results plotted in the latter form. Quadratic func- 
tions were found to represent the data satisfactorily. 

Since data was taken for 9 values of X and 7 values of 0, the two sets 
of functions were expressed as 

F^_ (0) = a.(X) + b^(X) 0 + c.(X). f A(12) 



where i = 1, 2 9 and a^, b^ and c^ are coefficients which depend on the 

value of X, and 

X.(F') = f (0) + g (0). [F’] + h (0).[F’l^ A(13) 

J o J J 9 J “ 

where j = 1, 2 7 and f.:, g. and h. are coefficients which depend on 

J J 

the value of 0. 

Next, the coefficients a^, b^ and c^ from the 9 equations represented 

by Eq. A (12) were approximated by a polynomial expansion in X, and the 

coefficients f ., g. and h. from the 7 equations represented by Eq. A(13) by 
J J J 



6 ? 



a polynomial expression in 0. Thus coefficients v/ere found for the expressions 



a. (X) = + A^X. + X^ + 

b. (X) = + B^X. + X^ + A(l4a) 



+ A(l4b) 



which gave the best fit to the 9 values of a^, b^, c^ and to the 7 values of 
f., g. and h.. This results in a single equation (Eq. A(12)) for FI (0) in 

t) J 0 r 

which the coefficients are analytic functions of X, and a single equation 

(Eq. A(15)) for X(F') in which the coefficients are analytic functions of 0. 

0 

The coefficients in Eq. A(l4) are given in Table A3- 

In application, the probe measurements give values to F^ and F 0 through 

Eq. A(10) and Eq. A(ll). Equation A(12) and Eq. A(13) must be solved for the 

unknown Ci and unknown X. An initial estimate of "the pitch angle was devised 
in the following way. A new function of the probe pressure differences was 
defined such that 



0. (X) = . c^x, . C 3 X. 



fj (0) = Fj, + + F 30 J 



g. (0) = . G^0j . 030" 



hj (0) - H3 . H^0j . H30^ 



F*(X,0) 

K 




A(15) 



When 0 was plotted as a function of F^ for each value of X, it was 

observed that the curves were nearly linear with similar slopes. Therefore 

the average of the 9 values of F' at each value of 0 was calculated and a 

K 

third order polynomial was found to fit the 7 average values when 0 was plotted 
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as a 



function the average F* (Fig. A9)* This gave 

\ 



0 (F') = J + K (F’) + L (F')^ f M (F')^ 

A. A A A 



A(16) 



where 0^ is the desired initial estimate of the pitch angle and F' is 
obtained from the probe measurements using Eq. A(15). The coefficients J, 

K, L & M are given in Table Ah. 

A2.U Application 

When applying the 5-hole probe to measure an unknown flow, the pressure 

differences read when p^ and p^ are balanced are reduced to the three factors 

F', FI and F' using Eq. A(10), Eq. A(11) and Eq. A(15). X and 0 are evaluated 
■ 0 7^ X 

using the calibration relations derived in the preceeding section. The yaw 
angle is read directly from the probe scale. The velocity is obtained from 
the product of X and which requires the input of the stagnation temperature. 

The reduction to velocity and pitch angle was programmed for the Hewlett- 
Packard Model 933 OA calculator. A flow chart of the program is shown in Fig. 
AlO. The program first calculates the initial estimate of 0 (=0^) using 
Eq. A(i 6). The initial estimate of X (=X^) is calculated from Eq. A(13) using 
0^ in Eq. A(lUb) to evaluate the coefficients f, g and h. The coefficients 
a, b and c in Eq. A(12) are then calculated using X^ in Eq. A(lUa). In 
general the value of F^ which could be calculated using these coefficients in 
Eq. A(12) would not be equal to the value calculated from the input data. An 

iterative procedure is begun using 0^ and X^ as initial values. 

In the adopted procedure, Eq. A(12) is viewed as an equation to be solved 

for an unknowi 0 in terms of a known F^ and known X. Equation A (13) is treated 

as an equation for X in terms of a kno\m 0. Equation A (12) is solved by fixed 
point iteration (9) beginning with the initial values 0^ and X^. In this 
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method, the (i+l)^^ value of 0 in the iteration is given in terras of the i^^ 
value by rewriting E:i. A (12) as 



- a - c 0 ^ 1/11 



A(17) 



For a fixed value of X(=X^), which fixes the coefficients a, b, and c, 
successive values of 0 are calculated \ising Eq. A(l7) until 
is satisfied, is a small quantity. A new value of X(=X^^^) is then 

calculated using 0 = 0^^^ ~ i'^ ^( 13 ) and the process repeated until 

|X^_i_j^ ~ ^i! *^x’ ®x ^ sraall quantity. 

The steps in the solution may be summarized as follows: 



1) Calculate 0 from Eq. A(i6) 

2) Calculate X^ from Eq. A (13) 

3) Calculate 0^^^.^^ from Eq. A(l?) 

4) If 10.. T - 0.1 <6/ do step 5)5 otherwise go to step 3) 

X* J- 1 p 

5 ) Calculate E'l* A(13) using 0 = 0^^^ = 0. 

6) If " ^i^ otherwise go to step 3)* 
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A . 3 Corrections to the Calibratio n or Application in an Annulus 
A3-1 Introduction 

The calibration described above was carried out where the presence of 
flow boundaries wouli not be expected to influenc*' the probe readings. A 
change in the calibration of the probe is expected when measurements are 
made within about ten probe diameters of a boundary, and when the probe size 
is significant compared to the flow passage area. A precise correction which 
would allow the use of the probe in any flow field is a difficult task and 
not attempted here. Instead, a calibration of the probe in an axial annular 
flow of similar geometry to the compressor annulus was carried out, and 
corrections to the "wall-free" calibration of section A2 were derived which 
were represented as a function of radial position. The test of the success 
of this proceedure was to compare measurements made in the compressor using 
two different types of probes. Agreement in the velocity profiles derived 
from the two probes would be expected if the corrections applied to each 
calibration were valid. Also, the distribution of true static pressure de- 
rived from the probe calibrations should be consistent with wall tap measure- 
ments on the tip and hub surfaces. In the following sections, the calibration 
tests are described and the derivation of corrections to the "wall-free" 
calibration is presented. The application of the revised calibration to 
compressor measurements is then given. 

A3 .2 Ca lib ration Test s in an Annulus 
An annular axial flow was generated as shown in Fig. All. An axisyra , 
metric body with a maximum diameter of 6"over 28 "of its length was supported 
in the center of a 10" pipe. At the upstream end, four thin airfoil struts 
were positioned ahead of the contoured nose section which generated a smooth 
contraction of tne Tlow passage, to a 2" annulus. 

r 
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Heavy struts supported the 



centerbody at the outlet flange where the flow exhausted into the laboratory. 
A view of the apparatus is shown in Fig. A12. Six equally spaced wall 
pressure taps were provided in the 10" jjipe at each of 3 axial stations. 
Corresponding pressure taps were provided on the centerbody surface. Pro- 
vision for radial probe surveys was ;tiade at the same three axial stations. 

It was established in preliminary surveys that the wall pressures at 
station 1 in Fig. All were uniform around the periphery to + 0.3^‘Ja of the 
average dynamic head. The outer (tip) wall pressures were consistently 
lower than inner (hub) wall pressures by 0.45^ of the dynamic head. The 
flow was measured to be axial to within the sensitivity of the 5-hole 
probe yaw angle measm'ement (+ .1°). Probe calibration was therefore carried 
out at station 1. Static pressure was assumed to vary linearly between hub 
and tip values, and the pitch angle was taken to be zero. 

Radial surveys were made from hub to tip at fixed flow rates using 
the 5-hole probe, a Kiel probe and a pitot-static probe in succession. 
Readings were taken at the same radial locations with each probe. Intervals 
of 0.2", decreasing to 0.05" near the walls, were used. Dynamic pressure was 
varied from 27 to 73 inches of water in 7-8 inch increments. The data for 
seven surveys is given in Table A2 . 

A3 -3 Analysis of the Results 

An example of the distribution of pressures measured by the three probes 
is given in Fig. A13. It was observed that the probes gave the same values 
of impact pressure at corresponding points to within the accuracy of the 
measurements except close to the outer wall where the velocity gradients were 
large. The static px'essure indicated by the pitot static probe was observed 
to vary very little across the annulus, and was in agreement with a linear 
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variation of the true static pressure between hub and tip measurements 
at the wall. Also, the wall pressures at hub and tip away from the traverse 

station were not measurably affected by the prove movement. Consequently, the 
5-hole probe indicated total presure (P]_) and a static pressure obtained by a 

linear interpolation between measurements at the walls were used to calculate 
the actual velocity (X) at the probe location. By applying the methods of 
section A2, the 5-hole probe pressure measurements were reduced to "wall-free" 
values of velocity and pitcri angle, (X and 0), which in general were different 
from the actual values, (X and0) at the probe locabiin. 

A velocity correction factor, K , was defined and examined as a function 
of position in the annulus and the indicated velocity, thus: 

K (H,X) = A (18) 

^ X 

whe're H = y/h, y is the distance outwards from the hub wall and h is the 
annulus depth. At each value, H^, of the probe position the dependence of 
on X was approxiiaated by a polynomial 

iC =q. +r.X+s.X^+ A (19) 

1 1 

where i = 1, 2, 12, and H. = 0.125, 0.15, 0.975. The coefficients 

in Eq. A (19) are listed in Table A5* 

The pitch angle at each probe position, H^, was found to be independent 
of X. A pitch correction factor, was defined which was then a function 
of position only: thus 

where p is a constant angle introduced to avoid division by zero. Since 
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A(21) 



0 = 0 in the calibration test, Eq. A(20) with = l80° becomes 

= - lo 



The values of K 



0 



approximated by a 



obtained for the tests in 
6^'"^ order polynomial, 



the annulus using Eq. A(2l) were 



(H) = H + z^ H^+ h3+z^h'^+ z^H^+ z^H^ A(22) 

to -995 root mean square. The coefficients in Eq. A(22) are listed in 
Table A6. 



A 3 . 4 App li c at i on 

It is assumed that the corrections required to the calibration free of 
wall effects when the 5-iiole probe is used to determine the flow field in a 
compressor annulus, are those established for an axial flow in the annulus 
of similar geometry. The probe readings are first reduced to vabes of X 
and 0 using the methods given in Section A2. Then, from Eq. .*.(l8) and 
Eq. A(20) with P = l80°, the corrected velocity and pitch angle are given by 



X = ^ ^(23) 

i - Kx 

and 

0 - (7^)- 0 ^ 180 *<84) 

At any probe position, H, the value of in Eq. A(24) is given by 
Eq. A(22) for which the coefficients are given in Table 6. 

The value of in Eq. A(23) must be obtained by interpolation. The 
uncorrected velocity X is kno\im at the probe position H; however, the correc- 
tion factor can be evaluated using Eq. A(19) only at discreet locations 
H^. The method adopted is a form of quadratic interpolation for uneven 
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intervals (10). First the value of i is found such that H is located in the 

interval hetween H. and H , as shown in Fig. A(14). The values of K , 

1 i+1 X..1 

K , K and K are calculated from Eq. A(19) using the coefficients 
^i ^i+1 ^i+2 

given in Table A5. The coefficients of quadratics for K passing through 
each of the two sets of three points H. , , H. , and H. , , H. 

1 - 1 ’ i ’ 1+1 1 1 + 1 ’ 1+2 

are calculated and each quadratic is used to calculate a value for at H. 
The arithmetic mean of the two values is taken for the interpolation. In 
the first or the last interval, or to extrapolate outside the given range, 
only one quadratic is used. 

The data reduction with corrections was programmed for the Hewlett- 
Packard Model 983OA calculator. The flow chart is given in Fig. A(l5). 

The steps are summarized as follows: 

1) Input probe pressure data, probe position, flow stagnation 
temperature, atmospheric pressure. 

2) Calculate ^ and X using the steps given in Section A2.4. 

3) Calculate with Eq. A (22). 

4) Calc\ilate K using quadratic interpolation between K values 

^i 

given by Eq. A(19). 

5) Calculate X and 0 using Eq. A(23) and Eq. A(24). 

This proceedure was applied to the calibration data as described in 
Section A3. 2 and listed in Table A2 . The results for one survey are shown 
for illustration in Table A?. At each probe position in the annulus, H, the 
uncorrected indicated velocity, X(UHMOD), the measured actual velocity, 
X(ACTUAL), and the corrected velocity, X(MOD) are compared. The largest 
error between the corrected and actual values was less than 0 . 3 % for all 
points in all surveys taken. 
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FIGURE A2, VIEV>f OF THE 5-HOLE PROBE SHOWING YAW 

ANGLE SCALE 
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FIGURE A3. SCHEMATIC OF THE APPARATUS FOR PROBE 
CALIBRATION FREE OF V/ALL EFFECT 




FIGURE A4. PITCH ANGLE BLOCK 
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FIGURE RE. VELDCITY PRDF 1 LE IN THE 

MER5URENENT PLRNE. 
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FIGURE RG. F»(THETR)Y5. PHI DHTR PDINT5 RNI> PDLYNDMIRL CURVE FITS. 
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FIGURE H7. F‘(PHI3 V5. PHI DATA POINTS RND POLYNUMIRL CURVE FITS. 
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figure: RB. X Y5 F’CTHFTR) DRTR BUNTS RND PDLYNDMIRL CURVE FITS. 
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FIGURE R3. PHI V5. RVERREE F'(LRMBDR) RND PDLYNDMIRL CURVE FITS. 
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FIGURE fll0. FLOW CHART FDR THE REDUCTION DF S-HDLE PROBE 
MEASUREMENTS TO VELOCITY AND PITCH ANGLE. 
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FIGURE All. SCHEMATIC OF THE APPARATUS FIR PRC3E 

CALIBRATION IN AN ANNULUS 




FIGURE A12. VIEW OF THE CALIBRATION APPARATUS 
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METHOD 




FIGURE AI5, FLOW CHART TO CORRECT PROBE CALIBRATION 



FOR MEASUREMENTS IN AN ANNULUS 
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TABLE A 1 .(C 0 NT 
SURV, BAROM. 
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TABLE A2. RESULTS OF PROBE CALIBRATION IN AN ANNULUS 
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TABLE A3. COEFFICIENTS FOR EQ. A(14) 
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TABLE A5. COEFFICIENTS FOR EQ. A (19) 
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TABLE A6. coefficients FOR EQ. A (22) 
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APPENDIX Bt DETERMINATION OF BLOCKAGE 



FACTOR FOR FLOW INTO THE ROTOR 
B.l. INTRODUCTION 

The analysis in Reference 1 established a 
blockage factor for the flow into the TRANSX com- 
pressor rotor which was used to correct inviscid 
mass flow calculations. The blockage factor at 
station 1 was defined as, 

Kg = ^m = measured flow rate B(l) 
inviscid flov/ rate 

and Kg was shown to be related to the displacement 
thickness of the boundary layer. Preliminary measure- 
ments were made ahead of the TRANSX compressor ro- 
tor to determine the boundary layer thickness at the 
outer wall [Reference l]. From the measurements, 
ignoring the hub wall boundary layer, and assuming 
that the inviscid flow field was described correct- 
ly, the result obtained was, 

= 0.989 

Then the flow into the rotor was described by the 

inviscid flow calculation for a total flow rate 

(iSr) which was related to the measured flow rate 

(\’z ) by the relation 
m 

* = • 

The purpose of the present study of the flow 
into the rotor was to measure the inviscid flow 
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• y. 














field, as well as the boundary layer region, and to 
relate the measurements to the analytical descrip- 
tion, It was necessary therefore to change the defin- 
ition of blockage factor, from that involving only 
boundary layers, to obtain the desired relationship. 

In the following sections, necessary parts of the 
analysis given in Reference 1 are repeated, and the 
method is given for relating the analytical distribu- 
tion of mass flux and flow angle to the distributions 
measured using the 5-hole probe, 

B.2. ANALYTICAL BACKGROUND AND METHOD 

Referring to Fig, Bl, and assuming that stream- 
lines are axisymmetric , continuity at station 1 
gives 

^ = J* fVcosX* (2irr) dr B(2) 

^i 

Non-dimensionalizing Eq, B(2), 

r 2 

W - ^ = r ° cosX d(£) B(3) 

^ rf 

The definition of the total flow function is given 

in Reference 7 as ^ 

$ = = X( l-x2) where 

ft t 

X = V/V. and V. = ^2c T. , 

t t P t 

V.’rjting R = r/^ , Eq , B(3) taken between R^ and R 2 
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becomes 



>2 



W = / J cos\ d(R'^) 

r.2 






Introducing the values of the flov/ function and flow 
angle at the outer wall, indicated by a "o" subscript, 
the above equation becomes 



W = S cos\^ / 
* o 





cos\ 

cd'sX 

o 



d(R^) 



In the TRANSX compressor ~ 0 , therefore 
2 

Rg _ ^ 

W = X 2 (“f”) B(4) 

^1 



The assumption is then made, in Reference 1 , 
that the pattern of streamlines does not depend on 
the flow rate, that ^ = ^/SqIs a function of R only 
and cos\ is another function of R. Based upon func- 
tions derived from a computer solution of the com- 
pressible flow, for Rj^ = R^ and R2 “ ^q* shown 

that 




cos\ d(R^) = 0.70836 



B( 4 a) 
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where, using Eq. B(4), 

. B(4b) 

Therefore the "theoretical" flow rate is proportional 
to the value of the analytical flow function at the 
outer wall. 

In order to relate the analytical description of 
the flov/ to the measured distributions of the flov/ 
function and angle, the "theoreticarl" mass flux was 
integrated and set equal to the integral of the 
measured distribution of the mass flux over that re- 
gion of the flow where the measurements were valid. 

The probe measurements were considered to be free of 
"immersion" or "wall" effects from 20 ^ to 80^ of the 
hub-to-tip survey. Therefore, using Eq, B(4) 



2 2 
r| 

S S cosX d(R^) = S I cosXd(R^) B(5) 

in 111 o Q 

R^ R^ 



defines a value of the analytical flow function at the 
outer wall v/hich corresponds to the measured flow 
conditions. Equation B(4b) then gave the "theoreti- 
cal" flov; rate corresponding to the measured flow rate, 
and Eq, B(l) was again used to evaluate a blockage 
factor. 
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B.3. PROCEDURE 



The analytical distribution of the total flow 
function [Ref. l] is represented as 

d-R'") = aj(l-§) >mj(l-$). B(6) 

where a^^, n^, are coefficients, listed in Table Bl, 
which were evaluated in Reference 1 from a computer 
solution using a curve fitting technique, f is un- 
knovm for a given R and must be determined as follov/s: 
1) Choose an initial value for § which does 
not exceed unity. Take 



2) 

3) 

4) 



= 1-.03(1-R^)^, R^ given. 

Set Y = l-l- , i=l, 2,..., n 
^ n. 

Calculate Rt = 1-a^Y "^-m^^Y 

2 

Test for convergence of calculated R ; 



2 2 
R^ - Rf 



5 ) 



<C t £= convergence limit. 
If yes, the solution is complete. 

If no, continue to step 5* 

Calculate the slope 



n^-1 



d(Rf)/dy = -a^n^Y -m^, and calculate 



j^i+1 - 



~ 2 2 
Rf-RT 
1 1 



d(Rf )/dy 



6) Repeat steps (2) through (5) until conver- 
gence occurs. 

Using this procedure the value of 1 could be calculated 
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for each probe position R 



A quantity was defined such that 




^ cosA.d(R^) 






B(7) 



was determined for a given and R 2 using numer- 
ical integration applied over unequal intervals. 

R^ was chosen to be the first point greater than 20^ 
of the distance from hub to tip and R 2 was chosen 
to be the last point before 80^ of the distance from 
hub to tip. 

The analytical distribution for the flow angle 
(X. ) used in Sq, E(4a) v.-as evaluated in Reference 1 



cosX = a 2 (1-R^) +m 2 (1-R^) B(8) 

where the coefficients are given in Table Bl, 

Using Eq, B(5), Eq, B(4b) and Eq, B(3) v/ith 
the perfect gas equation of state in Eq, B(7), 

■"o = ro — ■ J’ Pt '*('*■) ®(9) 

in G o 

Pt» and are given from measurements at each 
point in the probe survey and all other quantities in 
Eq, B(9) are constant for the survey. Since slight 
variations in the flow rate were noted during the sur- 
vey, a correction was introduced. The measured flow 
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rate, in Eq, B(9) was the average of data reduced 

from scans of fixed instrumentation taken before and 
after each survey. The value of corresponded to 
a differential pressure, Ah^, across a flow nozzle. 
Since the flow rate at any time was proportional to 
the square root of this differential pressure. Ah, 
the value of Ah was recorded at each survey point. 
Rearranging Eq, B(9)» with the correction introduced 



into the integration, the blockage factor is given 

t>y r 



Eq, B(6), B(7) and B(10) were programmed to calcu- 
late the blockage factor for flov/ at the rotor inlet 
from probe survey measurements. 

The correction for slight variations in the flow 
rate was also included in comparing the measured and 
analytical distributions of the mass flux (or flow 
function). The left hand side of Eq, B(5) is written, 
with the addition of the flow correction, as 




p*- 

“1 



B(10) 




^2 r^. 



B(ll) 



K 



m 



so that with Eq, B(7)» Eq» B(5) becomes 

K = $ K , 
m ^o 0 
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The analytical flow function, which was included 
on plots of the measured flow function, was obtained 
from 



$(R) = i (R) 



(8 



B(12) 



where J was obtained from Eq. B(6), was obtained 



from Eq, B(ll) and was already known from Eq, B(7). 
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FIGURE Bl. SCHEMATIC OF COMPRESSOR GEOMETRY 



TABLE Bl. CURVE FIT COEFFICIENTS 
FOR EQ. B(6) AND EQ. B(8) 



^EQUATION 


i 


a . 

1 


n. 

1 


m. 

1 


B(6) 


1 


1.79864 


0.38582 


22.3048 


B(8) 


2 


0.98927 


3.81679 


0.47149 
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